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Recent progress in the understanding of how externally driven magnetic reconnec¬ 
tion evolves is organized in terms of parameter space diagrams. These diagrams are 
constructed using four pivotal dimensionless parameters: the Lundquist number S', the 
magnetic Prandtl number Pm, the amplitude of the boundary perturbation iFo, and the 
perturbation wave number k. This new representation highlights the parameters regions 
of a given system in which the magnetic reconnection process is expected to be distin¬ 
guished by a specific evolution. Contrary to previously proposed phase diagrams, the 
diagrams introduced here take into account the dynamical evolution of the reconnection 
process and are able to predict slow or fast reconnection regimes for the same values 
of S and Pm, depending on the parameters that characterize the external drive, never 
considered so far. These features are important to understand the onset and evolution of 
magnetic reconnection in diverse physical systems. 


1. Introduction 


Magnetic reconnection is a process whereby the magnetic field line connectivity (New¬ 


comb 1958 Pegoraro 2012 Asenjo & Comisso 20151 is modified due to the presence 


of a localized diffusion region. This gives rise to a change in magnetic field line topol¬ 
ogy and a release of magnetic energy into kinetic and thermal energy. Reconnection of 
magnetic field lines is ubiquitous in laboratory, space and astrophysical plasmas, where 
it is believed to play a key role in many of the most striking and energetic phenomena. 


The most notable examples of such phenomena include sawtooth crashes (Yamada et al. 


1994 


1978 


Nicolas et al. 


20121 and major disruptions in tokamak experiments (Waddell et al. 


[Boozer 20121, solar and stellar flares (Masuda et al. 1994 Su et al. 20131, coro¬ 


nal mass ejections ( Lin fc Forbes [2000 Murphy et aL||20l2 i, magnetospheric substorms 
(0ieroset et al. 1|MI[ [Eastwood et a/.||2007), coronal heating ([Priest et al.||199^ |Cassai^ 


et al. 20081, and high-energy emissions in pulsar wind nebulae, gamma-ray bursts and 


jets from active galactic nuclei (Kagan et al. 2015 Sironi et al. 2015 Guo et al. 2015). 


An exhaustive understanding of how magnetic reconnection proceeds in various regimes 
is therefore essential to shed light on these phenomena. 

In recent years, for the purpose of organizing the current knowledge of the reconnection 
dynamics that is expected in a system with given plasma parameters, a particular form of 


phase diagrams have been developed (Ji & Daughton |2011 

Huang et a/.|2011 Daughton 

& Roytershteyn |20121|Huang & Bhattacharjee |2013[|Uassa 

c & Drake |20131 Karimabadi 

& Lazarian 2013 

). These diagrams classify what “phase” of magnetic reconnection should 
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occur in a particular system, which is identified by two dimensionless plasma parameters, 
the Lundquist number 


_ Ls'^A,u 

’-Dp ' 

(1.1) 

> 

III 

(1.2) 


and the macroscopic system size 


Here, indicates the system size in the direction of the reconnecting current sheet, va,u 
is the Alfven speed based on the reconnecting component of the magnetic field upstream 
of the diffusion region, = rjc^/Air is the magnetic diffusivity, and is the relevant 
kinetic length scale. This length scale corresponds to (see, e.g., Simakov fc Chacon [2008 
Comisso et al.|[2M3 ) 


Ik = 


di — cj UJpi 
Pr — Cg / 


for antiparallel reconnection, 
for guide-field reconnection. 


(1.3) 


Of course, ujpt is the ion plasma frequency, uJci is the ion cyclotron frequency, and Cg is 
the sound speed based on both the electron and ion temperatures. 


All the proposed phase diagrams |ji & Daughton 

2011 

Huang et aZ.|2011 

Daughton & 

Roytershteyn ||2012| Huang & Bhattacharjee ||2013 

Cassak & Drake ||2013 

Karimabadi 


& Lazarian 2013) exhibit a strong similarity and only a few minor differences. They 


are useful to summarize some of the current knowledge of the magnetic reconnection 
dynamics, but they lack fundamental aspects that can greatly affect the reconnection 


process (some caveats in the use of these diagrams have been discussed by Cassak fc Drake 
[ (20131). For example, they do not take into account the dependence of the reconnection 
process on the external drive or on the magnetic free energy available in the system. 
An attempt to include these effects has been discussed by Ji & Daughton (2011), who 


proposed to incorporate them by adjusting the definition of the Lundquist number, Eq. 
(O, but this solution should be viewed only as a rough way to circumnavigate the 
problem. A further issue is that these diagrams do not consider the evolution of the 
reconnection process and predict reconnection rates wich are always fast (the estimated 
reconnection inflow is always a significant fraction of va,u)- This, however, in not what 
is commonly observed in laboratory, space, and astrophysical plasmas, where magnetic 
reconnection exhibits disparate time scales and is often characterized by an impulsive 
behaviour, i.e., a sudden increase in the time derivative of the reconnection rate (see, 
e.g., Bhattacharjee ||2004 Yamada ]|2011 ). 

Here we propose a different point of view in which we include explicitly the effects 
of the external drive and the plasma viscosity (neglected in all previous diagrams) on 
the magnetic reconnection process by considering a four-dimensional parameter space. 
Then, in this four-dimensional diagram we identify specific domains of parameters where 
the reconnection process exhibits distinct dynamical evolutions. In other words, in each 
of these domains the reconnection process goes through diverse phases characterized 
by different reconnection rates. This analysis leads us to evaluate in greater detail the 
dynamical evolution of a forced magnetic reconnection process, while collisionless effects 
have not been taken into account in the present work. We introduce the considered model 
of forced magnetic reconnection in Sec. whereas Sec. is devoted to the presentation 
of the possible evolutions of the system and the conditions under which these different 
evolutions occur. In Sec. we construct the parameter space diagrams that show which 
reconnection evolution is expected in a system with given characteristic parameters. 
Finally, open issues are discussed in Sec. 
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Figure 1. Geometry of the Taylor model. The equilibrium magnetic field component By is 
sheared in the x direction, being null at x = 0. The plasma is bounded by perfectly conducting 
walls at X = ±L, while it is periodic in the y direction. Magnetic reconnection is driven at x = 0 
by the perturbation Ho cos{ky) at the perfectly conducting walls. 


2. Forced magnetic reconnection in Taylor’s model 

Magnetic reconnection in a given system is conventionally categorized as spontaneous 
or forced. Spontaneous magnetic reconnection refers to the case in which the reconnection 
arises by some internal instability of the system or loss of equilibrium, with the most 
typical example being the tearing mode. Forced magnetic reconnection instead refers to 
the cases in which the reconnection is driven by some externally imposed flow or magnetic 
perturbation. In this case, one of the most important paradigms is the so-called “Taylor 
problem”, which consists in the study of the evolution of the magnetic reconnection 
process in a tearing-stable slab plasma equilibrium which is subject to a small amplitude 
boundary perturbation. This situation is depicted in Fig.j^ where the shared equilibrium 
magnetic field has the form 

B = B^e^ + {x/L)Boey, (2.1) 

with Bz, Bq and L as constants, and the perfectly conducting walls which bound the 
plasma are located at x = ±L. Magnetic reconnection is driven at the resonant surface 
X = 0 by a deformation of the conducting walls such that 


±L =F -0 cos{ky ), 


( 2 . 2 ) 


where k = ‘I'k fLy is the perturbation wave number and Hg is a small (<?; L) displacement 
amplitude. The boundary perturbation is assumed to be set up in a time scale that is 
long compared to the Alfven time ta = L/va, with va = B^jy/A t: but short compared 
to any characteristic reconnection time scale. Hence, the plasma can be considered in 
magnetostatic equilibrium everywhere except near the resonant surface at x = 0. 

The first and probably most important contribution to unveiling the behaviour of 


forced magnetic reconnection in Taylor’s model is due to Hahm & Kulsrud (1985), who 


showed that very small amplitude boundary perturbations cause an initial linear phase 
in which a current sheet builds up at the resonant surface, and successive phases in which 
the reconnection process evolves according to a linear resistive regime and a nonlinear 


Rutherford regime (Rutherford 1973). The scenario discussed by Hahm & Kulsrud 


(1985), which is characterized by a very slow evolution of the reconnection process, was 
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complemented some years later by Wang & Bhattacharjee (19921, who showed that 


larger perturbations may foster reconnection to proceed through the nonlinear regime 


according to a Sweet-Parker-like evolution (Waelbroeck 19891, which only on the long 


time scale of resistive diffusion gives way to a Rutherford evolution. The scenario out¬ 


lined by Wang & Bhattacharjee (1992) is characterized by a reconnection evolution 


faster than that presented by Hahm & Kulsrud (19851, but it could still be slow for very 
small values of plasma resistivity, since in both the Sweet-Parker-like (Waelbroeck [1989 ) 


and Rutherford (Rutherford |1973) regimes, the reconnection rate is strongly dependent 


on the resistivity, which is known to be extremely small in many laboratory fusion plas¬ 
mas and space/astrophysical plasmas. However, recent works (Comisso et aL|2014 2015) 


have shown that relatively large boundary perturbations lead to a different reconnection 
evolution in plasmas with small resistivity and viscosity. In these cases, after a linear in¬ 
ertial phase and an initial nonlinear regime characterized by a gradually evolving current 
sheet, the reconnection suddenly enters into a fast reconnection regime distinguished by 
the disruption of the current sheet due to the development of secondary magnetic islands 
(usually called plasmoids ( Biskamp]|2000 Loureiro et al.|[2007|). 


In addition to the works discussed above, which adopt a magnetohydrodynamic (MHD) 
description of the plasma, we emphasize that many other efforts have been devoted to 
investigate the Tayor problem assuming MHD, two-fluid and kinetic descriptions (see 


Wang & Bhattacharjee 1992| Ma et al 1996 Rem & Schep 1998| Vekstein & Jain | 

1998 Avinash et aL||1998 Vekstein & Jain 

1999 Valori et al.|20001 Fitzpatrick ||2003 

Fitzpatrick et al.\ 

20031 Fitzpatrick ||2004a|b[ 

Cole & Fitzpatrick ||2004 

|Bian & Veksteinj 

20051 |Birn et al. 

|2005t IVekstein & Bian 2006| |Birn & Hesse ||2007[ | 

i’itzpatrick ||2008t 


Hosseinpour fc Vekstein 2008[ 


Hosseinpour [2013 |Dewar et al. 


Gordovskyy et al. 2010a|b Lazzaro & Comisso ]|0TTt 


2013). Indeed, the Taylor problem has important appli¬ 


cations besides being interesting from the point of view of basic physics. For instance, 
in laboratory fusion plasmas the Taylor model represents a convenient way to study 
magnetic reconnection processes driven by resonant magnetic perturbations, while in as- 
trophysical plasmas this model can be adopted to study magnetic reconnection forced by 
the motions of photospheric flux tubes. 


3. Evolution of the reconnection process in Taylor’s model 

In this section we review the present understanding of the forced magnetic reconnection 
dynamics in Taylor’s model focusing on a visco-resistive plasma with greater than 1. 
As shown by Hahm & Kulsrud (1985), this dynamics always starts with a linear inertial 
phase in which a current sheet builds up at the resonant surface and shrinks inversely 
in time. Concurrently, the current density at the Ai-point increases linearly in time. The 
reconnection rate during this phase can be evaluated by recalling that the current density 
is proportional to the out-of-plane electric field at the Ai-point, which is equal to 




t 


0 ^ 0 " 




(3.1) 


for t (TAjkL'f'l'^ (Fitzpatrick 


2003 


Comisso et al. 


2015). Here, i/’ stands for 


the magnetic flux function of the perturbed magnetic field in the reconnection plane 
{5B± = \/ip X e^), Ti, = Li^/v and indicate the characteristic time for 

viscous and resistive diffusion, respectively, while A(, = 2fc/sinh(fcL) parametrizes the 
contribution of the external source perturbation to the gradient discontinuity of the 
magnetic flux function at the resonant surface. It is important to point out that this 
phase is characterized by a non-constant-^ behaviour of the magnetic flux function across 
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the island. Depending on whether or not this property persists until the beginning of the 
nonlinear regime, different scenarios may occur. 


3.1. Hahm-Kulsrud scenario 

If the boundary perturbation is such that ( Fitzpatrick [2003 Comisso et aL|[2M5 | 

KkL) At “ ' 


5'n = B, 


0^0 


«(- 


U ‘7)) 




(3.2) 


after the inertial phase the reconnection process evolves trough a visco-resistive phase, 
which is a linear regime characterized by a constant-'^ behaviour, i.e., the perturbed 
magnetic flux function can be treated as a constant in x over the width of the reconnection 
layer. During this phase the reconnection rate is given by 

(3.3) 

where Ag = 2fc/tanh(fcL) is the standard tearing stability parameter and r* is a charac¬ 
teristic time defined as ( Fitzpatrick ||2003] [Comisso et al.||2015 ) 

5/6 


= ^g2/3 ^ (e) 


1/3 


r(i) 


1/6 KkLJ 


(3.4) 


with F indicating the Gamma function. Eq. (3.3) is valid for t {TAjkL'f'l^ 

(Fitzpatrick [2003 Comisso et al.|[20T5l and a magnetic island width much smaller than 


the linear layer width, i.e., w (Porcelli 


1987 


Fitzpatrick 


1993). If the perturbation is sufficient to drive the magnetic island into the nonlinear 
regime {w > ^i/tj), the visco-resistive phase ends up into a Rutherford evolution, whose 
island width growth is governed by the Rutherford equation 

^'o 


I 


Bydw _ , 

L2 dt ~ ° 


A' 




(3.5) 


where X = 0.823 ( Rutherford [1973 Fitzpatrick [1993 ). This is a very slow reconnection 
evolution in which the reconnection rate can be evaluated analytically in the two limits 
( Hahm fc Kulsrucl~||1985 Comisso et aL|[20T5 ) 

/ 3t 


dtfjlx = 


2A(^'n 


{—^o)tnl \tnlJ 


for t < tnl , 


dt'iplx = 


2A;^o 


tanh 


(-) 

X'rNLj 


cosh 


(M 

\tnlJ 


for t > TxL 


where 


41 


tnl = 


a: Sr 


(-a')lV(-a') l 


1/2 


(3.6) 


(3.7) 


(3.8) 


3.2. Wang-Bhattacharjee scenario 

If the boundary perturbation is such that ( Fitzpatrick [2003 Comisso et aL|[2M5 | 


(3.9) 


the non-constant-i/) behaviour characteristic of the inertial phase lingers until the nonlin¬ 
ear regime is entered. Therefore, since in this case the magnetic island grows faster than 
the current can diffuse out of the reconnecting layer, the evolution of the reconnection 
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process is distinguished by a strong current sheet at the resonant surface (Waelbroeck 


1989). The reconnecting current sheet turns out to be stable if the boundary perturbation 
is such that ( Comisso et aL|[20T5 l 

1/2 


^0 — Bq^q < CBqL 


k ta 


1 + 


= 'fr 


(3.10) 


where the multiplicative constant C ~ 2e“^ depends on the critical inverse aspect ratio 
of the reconnecting current sheet (specified later). In this case, the reconnection process 
follows a Sweet-Parker evolution (modified by plasma viscosity ( Park et a/.|1984 |), whose 
reconnection rate in Taylor’s model is ( Comisso et a/.|[2015 ) 


dt^\ 


X 


^BoLiX 




1 + ^ 
Tu 


(3.11) 


Finally, the Sweet-Parker type of evolution gives way to a Rutherford evolution on the 
time scale of resistive diffusion. 


3.3. Our scenario 

If the boundary perturbation satisfies ( Fitzpatrick ||2003 Comisso et a/.|[20T5 ) 


'W 


and also the condition (Comisso et a/.|[2015) 


^o>CBoL-^:^ 

^ c. In 


1 + Bl 

T,, 


1/2 




(3.12) 


(3.13) 


the reconnection process does not reach a stable Sweet-Parker regime, but a different 
situation occurs. A gradually thinning current sheet evolves until its aspect ratio reaches 
the limit that allows the plasmoid instability to develop. The growth of the plasmoids 
leads to the disruption of the current sheet, and therefore to a dramatic increase of the 
reconnection rate. The reconnection rate during this plasmoid-dominated phase has been 
evaluated in a statistical steady state as ( Comisso et aZ.|[20T5 ) 

- 1/2 


'■3' l2_-l 


dtipp « ecBoL{A^^o) T 


1 + A- 

Tu 


(3.14) 


where Cc = dcjLc is the critical inverse aspect ratio of the reconnecting current sheet. 
This quantity, whose value has been found to lie in the range 1/100 - 1/200 by means 


of numerical simulations (|Bhattacharjee et a/.||2009 

Samtaney et aL||2009t Cassak et al. 

to 

o 

o 

Skender & Lapenta |2010 

Huang & Bhattacharjee 120101, represents the threshold 


below which the reconnecting current sheet becomes unstable to the plasmoid instability 


( Loureiro et aL|[2007 ). 


4. Phase diagrams 

In this section we illustrate the domain of existence of the three different scenarios 
described before with the help of appropriated parameter space maps. For the sake of 
clarity we state again the three type of reconnection evolutions we are referring to: 


(1) Hahm-Kulsrud scenario ( Hahm &: Kulsrud~||1985 Fitzpatrick [2003 ), 
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(2) Wang-Bhattacharjee scenario ( Wang &: Bhattacharjee |[T992 Fitzpatrick ||2003 |, 

(3) Our scenario (Comisso et a7|[20T5|. 


Since each of these scenarios includes different phases/regimes of reconnection, the con¬ 
cept of “phase diagrams” is intended here in a broader sense. Due to this fact, they could 
also be defined in a more general way as “scenario diagrams”. This kind of diagrams can 
be constructed from the conditions summarized in the previous section. Therefore, the 
possible evolutions of the reconnection process may be organized in a four-dimensional 
parameter space map with t]>g = To/SgT, k = kL, S = LvA/Drj, and = v on the 
four axes. However, due to the difficulty in the visualization of such a four-dimensional 
diagram, it is convenient to consider two-dimensional slices for fixed values of two of the 
four parameters. 

Let us first consider four two-dimensional slices with fixed values of the magnetic 
Prandtl number and perturbation wave number. Assuming that the Hahm-Kulsrud sce¬ 
nario (which occurs if dto ’^w) may hold until d>o = 4'w/3, the corresponding dia¬ 
grams for (a) it = 1/8, Pm = 5, (b) k = 1/8, Pm = 500, (c) k = 2, Pm = 5, (d) fc = 2, 
Pm = 500, are shown in Figs. da) - gd). From this plots it is clear that the Wang- 
Bhattacharjee scenario is limited to a small range of values of the Lundquist number 
and the source perturbation amplitude. Increasing values of the the magnetic Prandtl 
number and perturbation wave number extend the domain of existence of this possi¬ 
ble type of evolution of the system. However, after a threshold value of the Lundquist 
number (identified by the intersection of the two black lines representing = ^w/3 
and 'ho = ^c)) the Wang-Bhattacharjee scenario cannot occur because it is not possible 
to obtain a stable Sweet-Parker-type evolution. In these cases, the Hahm-Kulsrud sce¬ 
nario is facilitated by very small perturbation amplitudes, whereas larger perturbations 
lead the system to a fast reconnection regime as described in Sec. s Note that while 


previously proposed phase diagrams always predict fast reconnection (Ji & Daughton 

2011 Huang et a/.|2011 

Daughton & Roytershteyn |2012| Huang & Bhattacharjee |2013 

Cassak & Drake |2013 

<arimabadi & Lazarian |2013 

1, in clear contrast to what happens 


in nature, our diagrams show that reconnection proceeds very slowly (region (1)) if the 
source perturbation is not sufficiently large. 


Let us now examine the effect of the plasma viscosity by considering the domain of 
existence of the different scenarios as a function of the parameters and Pm- Fig.j^a) 
shows the functions 'I'o = 'I'rv/3 and 'Fg = for S = 10® and k — 0.5. For 'I'c < 
the threshold for the plasmoids formation coincides with that for the nonlinear evolution 
characterized by a strong reconnecting current sheet. Therefore, for < 'Fw/3 an 
increase of the amplitude perturbation dtg drives the system directly from scenario (1) 
to scenario (3). This situation is depicted in Fig. [^b), where it is clearly shown that the 
increase of the magnetic Prandtl number has the effect of making possible or extending 
the domain of existence of scenario (2). 

To clarify the effect of the plasma viscosity we also delineate the boundaries of the 
diverse evolutions (l)-(3) in a parameter space map (as in Fig. for fixed 

k = 0.5 but different values of Pm- This is shown in Fig. 0 for = 5 - 5 X 10®. The 
increase of the magnetic Prandtl number extends the domain of existence of the slow 
reconnection scenario (1) at the expense of the fast reconnection scenario (3). The area 
of existence of scenario (2) remains almost unchanged, but shifted towards higher values 
of the Lundquist number. 

We now examine in more detail how the possible evolutions of the forced magnetic 
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Figure 2. Two-dimensional slices of a phase/scenario diagram for forced magnetic reconnection 
in the magnetohydrodynamical Taylor model. Fixed parameters are (a) fc = 1/8, Pm = 5, (b) 
fc = 1/8, Pm = 500, (c) fc = 2, Pm = 5, and (d) fe = 2, Pm = 500. The numerical labels indicate 
the (1) Hahm-Kulsrud scenario, (2) Wang-Bhattacharjee scenario, and (3) our scenario. The 
boundaries between the different scenarios are identified by the functions = ^iu/3 and 
for 




(d) 




reconnection process depend on the wave number of the boundary perturbation. Fig. 
l^a) shows the thresholds 4'o = '^vu/3 and = 'I'c as a function of k for fixed values 
of S' = 10® and Pm = 5. Below a critical perturbation wave number fc* (corresponding 
to fc* « 1 for the fixed parameters used in Fig. iJa)), every time the non-constant-'0 
magnetic island pass into the nonlinear regime, the evolution of the system leads to the 
plasmoid-dominated phase predicted in scenario (3). The domains of existence of the 
possible evolutions (l)-(3) are illustrated in Fig. [^b). The scenario discussed by Wang 
and Bhattacharjee happens only for a small range of parameters. Not also that 

scenario (3) is facilitated for k < k*, while scenario (1) may occur for large amplitude 
boundary perturbations if fc ^ fc*. 
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Figure 3. (a) Thresholds Tvu/3 (red line) and 'he (blue line) as a function of the magnetic 
Prandtl number Pm. for S = 10®, k — 0.5 and C = 2(150)^. (b) Corresponding two-dimensional 
slice of the phase/scenario diagram identifying the (1) Hahm-Kulsrud scenario, (2) Wang-Bhat- 
tacharjee scenario, and (3) our scenario. 



Pm=50 

— P,„=500 

— P,„=5000 

— P,„=50000 


Figure 4. Boundaries (identified by the functions To = Tvy/3 and To = Tc for Tc > Tu//3) 
of the different possible evolutions of the reconnection process for k = 0.5 and various values of 
the magnetic Prandtl number. 


To better evaluate the effects of k on the possible evolutions of the reconnection process, 
we plot in Fig. the boundaries between the scenarios (l)-(3) in a parameter space map 
(To, S') (as in Figs. I and 1^ for fixed Pm = 5 but different values of k. The maximum 
area of existence of scenario (2) occurs for fc ~ 1, while for fc <C 1 and fc ^ 1 the scenario 
(2) appears for a very limited range of (Tq, k) parameters. Note also that the scenario (3) 
is greatly facilitated in the case of very large perturbation wave numbers (fc ^ 1) while 
the scenario (3) is facilitated by relatively large amplitude perturbations with {k < 1). 
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Figure 5. (a) Thresholds Tvu/3 (red line) and ^'c (blue line) as a function of the perturbation 
wave number k for S = 10®, Pm = 5 and C = 2(150)^. (b) Corresponding two-dimensional slice of 
the phase/scenario diagram identifying the (1) Hahm-Kulsrud scenario, (2) Wang-Bhattacharjee 
scenario, and (3) our scenario. 



— k=0.02 

— k=0.04 

— k=0.08 

— k=0.16 

— k=0.32 

— k=0.64 
k=1.28 
k=2.56 

— k=5.12 

— k=10.24 


Figure 6. Boundaries (identified by the functions 4'o = ^vu/3 and 'I'o = (fc for 4'c > 4'w/3) 
of the different possible evolutions of the reconnection process for Pm = 5 and various values of 
the perturbation wave number. 


5. Discussion 


The introduction of a new type of phase/scenario diagrams that include explicitly the 
effects of the external drive has allowed us to graphically organize in a detailed way 
the possible evolutions of forced magnetic reconnection processes in collisional plasmas. 


In contrast to previous versions of the phase diagrams (Ji & Daughton |2011 Huang 

\et at. 

201 1| Daughton & Roytershteyn 

2012| [Huang & Bhattacharjee 

20131 ICassak & 

|Drake 

2013 

Karimabadi & Lazarian 1120131, this new representation highlights regions 


of the parameter space ('I'oj^) ‘S', Pm) in which reconnection is a slow diffusive process 
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(Sec. 3.1) in addition to regions where reconnection can be fast (Secs. 3.2 and 3.3). We 
recall that by fast we mean that the out-of-plane inductive electric field at the X-point 
is a significant fraction of the one evaluated upstream of the reconnection layer. We 
also emphasize that this kind of diagrams respond to the criticism moved by [Cassak 
& Drake (2013) concerning the fact that these diagrams are not able of taking into 


account the dynamical evolution of the reconnection process from a slow to a fast regime 
inside of a given region of the parameter space. Indeed, scenarios (l)-(3) describe the 
forced magnetic reconnection process from the current sheet formation all the way to 
their specific nonlinear evolution. 

We would like to remark that while the proposed parameter space diagrams represent 
a valid way to summarize the current knowledge of the forced magnetic reconnection 
dynamics in a collisional plasma, there are a number of conditions that may significantly 
affect the reconnection process which have not been addressed in this paper. For instance, 
two-fluid/kinetic effects should be considered if the length scale associated with the width 
of the reconnecting current sheet becomes of the order or smaller than the characteristic 
length scales of these effects. In fact, in antiparallel reconnection (i.e., in the absence of a 
guide magnetic field). Hall effects (Birn et al.|200I Simakov fc Chacon [2008 ) are known 


to enhance the reconnection rate, as well as effects associated to finite electron inertia 


Ottaviani & Porcelli |1993 

Comisso & Asenjo 

2014 

), electron pressure ( 

Kleva et al. 


1995 Grasso et al. 1999) and ion gyration (Comisso et al. 2013) are known to increase 


the reconnection rate in the presence of a strong guide field. We would like to remark 
also that a common condition in many physical systems is the presence of velocity flows. 


which are known to suppress the reconnection 

Fitzpatrick |1993 

Waelbroeck et al. 2012 

or to alter the reconnection rate ( 

Cassak 2011 

Tassi et al. 2014) 

In this case our analysis 


should be extended by considering also the effects of a plasma flow on the reconnection 


dynamics. Similarly, also the effects of turbulence should be considered ( Servidio et al. 
2009 Karimabadi fc Lazarian ||2013 ) in order to obtain a more complete description of 
the magnetic reconnection dynamics. 

Finally, it is important to recall that all the presented diagrams of magnetic reconnec¬ 
tion are based on two-dimensional models and simulations. At present, the knowledge of 
how magnetic reconnection evolve in large three-dimensional systems is still far behind 
our understanding of what happens in two-dimensional systems. Therefore, despite the 
great progress achieved in recent years |Borgogno et a/.|[2005 Yin et al(]|2008 Daughton 
et al. ][ml| |Wyper fc Pontin |2014[ ), other work is needed in this direction before we can 
implement a phase diagram description of three-dimensional magnetic reconnection. 


The authors would like to acknowledge fruitful conversations with Richard Fitzpatrick 
and Enzo Lazzaro. This work was carried out under the Contract of Association Euratom- 
ENEA and was also supported by the U.S. Department of Energy under Contract No. 
DE-FG02-04ER-54742. 


REFERENCES 

Asenjo F. a. & Comisso L. 2015 Phys. Rev. Lett. 114, 115003. 

Avinash K., Bulanov S. V., Esirkepov T., Kaw P., Pegoraro F., Sasorov P. V. & Sen 
A. 1998 Phys. Plasmas 5, 2849. 

Bhattacharjee a. 2004 Annu. Rev. Astron. Astrophys. 42, 365. 

Bhattacharjee a., Huang Y.-M., Yang H. & Rogers B. 2009 Phys. Plasmas 16, 112102. 
Bian N. & Vekstein G. 2005 Phys. Plasmas 12, 072902. 

Birn J., Drake J. F., Shay M. A., Rogers B. N., Denton R. E., Hesse M., Kuznetsova 


































































12 


L. Comisso, D. Grasso and F.L. Waelbroeck 


M., Ma Z. W., Bhattacharjee A., Otto A. & Pritchett P. L. 2001 J. Geophys. Res. 
106, 3715. 

Birn J., Galsgaard K., Hesse M., Hoshino M., Huba J., Lapenta G., Pritchett P. L., 
Schindler K., Yin L., Buchner J., Neukirch T. & Priest E. R. 2005 Geophys. Res. 
Lett. 32, L06105. 

Birn J. & Hesse M. 2007 Phys. Plasmas 14, 082306. 

Biskamp D. 2000 Magnetic Reconnection in Plasmas, Cambridge University Press. 

Boozer A. H. 2012 Phys. Plasmas 19, 058101. 

Borgogno D., Grasso D., Porcelli F., Califano F., Pegoraro F. & Farina D. 2005 
Phys. Plasmas 12, 032309. 

Gassak P. a., Mullan D. J. & Shay M. A. 2008 Astrophys. J. 676, L69. 

Gassak P. a., Shay M. A. & Drake J. F. 2009 Phys. Plasmas 16, 120702. 

Gassak P. A. 2011 Phys. Plasmas 18, 072106. 

Gassak P. A. & Drake J. F. 2013 Phys. Plasmas 20, 061207. 

Gole a. & Fitzpatrick R. 2004 Phys. Plasmas 11, 3525. 

Gomisso L., Grasso D., Waelbroeck F. L. & Borgogno D. 2013 Phys. Plasmas 20, 092118. 
Gomisso L. & Asenjo F. A. 2014 Phys. Rev. Lett. 113, 045001. 

Gomisso L., Grasso D. & Waelbroeck F. L. 2014 J. Phys.: Conf. Ser. 561, 012004. 
Gomisso L., Grasso D. & Waelbroeck F. L. 2015 Phys. Plasmas 22, 042109. 

Daughton W., Roytershteyn V., Karimabadi H., Yin L., Albright B., Bergen B. & 
Bowers K. 2014 Nature Phys. 7, 539. 

Daughton W. & Roytershteyn V. 2012 Spaee Sci. Rev. 172, 271. 

Dewar R. L., Bhattacharjee A., Kulsrud R. M. & Wright A. M. 2013 Phys. Plasmas 
20, 082103. 

Eastwood J. P., Phan T., Mozer F. S., Shay M. A., Fujimoto M., Retino A., Hesse 
M., Balogh a., Lucek E. A. & Dandouras I. 2007 J. Geophys. Res. 112, 6235. 
Fitzpatrick R. 1993 Nucl. Fusion 33, 1049. 

Fitzpatrick R. 2003 Phys. Plasmas 10, 2304. 

Fitzpatrick R., Bhattacharjee A., Ma Z. W. & Linde T. 2003 Phys. Plasmas 10, 4284. 
Fitzpatrick R. 2004a Phys. Plasmas 11, 937. 

Fitzpatrick R. 2004b Phys. Plasmas 11, 3961. 

Fitzpatrick R. 2008 Phys. Plasmas 15, 024503. 

Gordovskyy M., Browning P. K. & Vekstein G. E. 2010a Astron. Astrophys. 519, A21. 
Gordovskyy M., Browning P. K. & Vekstein G. E. 2010b Astrophys. J. 720, 1603. 
Grasso D., Galifano F., Pegoraro F. & Porcelli F. 1999 Plasma Phys. Control. Fusion 
41, 1497. 

Guo F., Liu Y.-H., Daughton W. & Li H. 2015 submitted to Astrophys. J., arXiv:1504.02193 
[astro-ph.HE]. 

Hahm T. S. & Kulsrud R. M. 1985 Phys. Fluids 28, 2412. 

Hosseinpour M. & Vekstein G. 2008 Phys. Plasmas 15, 022904. 

Hosseinpour M. 2013 J. Plasma Physies 79, 519. 

Huang Y.-M. & Bhattacharjee A. 2010 Phys. Plasmas 17, 062104. 

Huang Y.-M., Bhattacharjee A. & Sullivan B.P. 2011 Phys. Plasmas 18, 072109. 
Huang Y.-M. & Bhattacharjee A. 2013 Phys. Plasmas 20, 055702. 

Ji H. & Daughton W. 2011 Phys. Plasmas 18, 111207. 

Kagan D., Sironi L., Gerutti B. & Giannios D. 2015 Space Sci. Rev. 186, 1. 

Karimabadi H. & Lazarian A. 2013 Phys. Plasmas 20, 112102. 

Kleva R. G., Drake J. F. & Waelbroeck F. L. 1995 Phys. Plasmas 2, 23. 

Lazzaro E. & Gomisso L. 2011 Plasma Phys. Control. Fusion 53, 054012. 

Lin j. & Forbes T. G. 2000 J. Geophys. Res. 105, 2375. 

Loureiro N. F., Schekochihin A. A. & Cowley S. C. 2007 Phys. Plasmas 14, 100703. 

Ma Z. W., Wang X. & Bhattacharjee A. 1996 Phys. Plasmas 3, 2427. 

Masuda S., Kosugi T., Hara H., Tsuneta S. & Ogawara Y. 1994 Nature 371, 495. 
Murphy N. A., Miralles M. P., Pope C. L., Raymond J. C., Winter H. D., Reeves K. 
K., Seaton D. B., van Ballegooijen A. A. & Lin J. 2012 Astrophys. J. 751, 56. 


Phase Diagrams of Forced Magnetic Reconnection in Taylor’s Model 


13 


Newcomb W. A. 1958 Ann. Phys. 3, 347. 

Nicolas T., Sabot R., Garbet X., Lutjens H., Luciani J.-F., Guimaraes-Filho Z., 
Decker J. & Merle A. 2012 Phys. Plasmas 19, 112305. 

0EROSET M., Phan T. D., Fujimoto M., Lin R. P. & Pepping R. P. 2001 Nature 412, 414. 
Ottaviani M. & Porcelli F. 1993 Phys. Rev. Lett. 71, 3802. 

Park W., Monticello D. A. & White R. B. 1984 Phys. Fluids 27, 137. 

Pegoraro F. 2012 Europhys. Lett. 99, 35001. 

Porcelli F. 1987 Phys. Fluids 30, 1734. 

Priest E. R., Foley G. R., Heyvaerts J., Arber T. D., Culhane J. L. & Acton L. W. 
1998 Nature 393, 545. 

Rem j. & SCHEP T. J. 1998 Plasma Phys. Controlled Fusion 40, 139. 

Rutherford P. H. 1973 Phys. Fluids 16, 1903. 

Samtaney R., Loureiro N. F., Uzdensky D. A., Schekochihin A. A. & Gowley S. G. 
2009 Phys. Rev. Lett. 103, 105004. 

Servidio S., Matthaeus W.H., Shay M.A., Gassak P.A. & Dmitruk P. 2009 Phys. Rev. 
Lett. 102, 115003. 

Skender M. & Lapenta G. 2010 Phys. Plasmas 17, 022905. 

SiMAKOV A. & Ghacon L. 2008 Phys. Rev. Lett. 101, 105003. 

SiRONi L., Petropoulou M. & Giannios D. 2015 Mon. Not. R. Astron. Soc. 450, 183. 

Su Y., Veronig a. M., Holman G. D., Dennis B. R., Wang T., Temmer M. & Gan W. 
2013 Nature Phys. 9, 489. 

Tassi E., Grasso D. & COMISSO L. 2014 Fur. Phys. J. D 68, 88. 

Valori G., Grasso D. & de Blank H. J. 2000 Phys. Plasmas 7, 178. 

Vekstein G. & Jain R. 1998 Phys. Plasmas 5, 1506. 

Vekstein G. & Jain R. 1999 Phys. Plasmas 6, 2897. 

Vekstein G. & Bian N. 2006 Phys. Plasmas 13, 122105. 

Waddell B. V., Garreras B., Hicks H. R., Holmes J. A. & Lee D. K. 1978 Phys. Rev. 
Lett. 41, 1386. 

Waelbroeck E. L. 1989 Phys. Fluids B 1, 2372. 

Waelbroeck F. L., Joseph L, Nardon E., Becoulet M. & Fitzpatrick R. 2012 Nucl. 
Fusion 52, 074004. 

Wang X. & Bhattacharjee A. 1992 Phys. Fluids B 4, 1795. 

Wang X. & Bhattacharjee A. 1992 Astrophys. J. 401, 371. 

Wyper P. F. & Pontin D. I. 2014 Physics of Plasmas 21, 082114. 

Yamada M., Levinton F. M., Pomphrey N., Budny R., Manickam J. & Nagayama Y. 

1994 Phys. Plasmas 1, 3269. 

Yamada M. 2011 Phys. Plasmas 18, 111212. 

Yin L., Daughton W., Karimabadi H., Albright B. J., Bowers K. J. & Margulies J. 
2008 Phys. Rev. Lett. 101, 125001. 



